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Abstract 

We study salt-induced charge overcompensation and charge inversion of flexible polyelectrolytes 
via computer simulations and demonstrate the importance of ion excluded volume. Reentrant 
condensation takes place when the ion size is comparable to monomer size, and happens in a 
middle region of salt concentration. In a high-salt region, ions can overcharge a chain near its 
surface and charge distribution around a chain displays an oscillatory behavior. Unambiguous 
evidence obtained by electrophoresis shows that charge inversion does not necessarily appear with 
overcharging and occurs when the ion size is not big. These findings suggest a disconnection of 
resolubilization of polyelectrolyte condensates at high salt concentration with charge inversion. 

PACS numbers: 82.35.Rs, 36.20.Ey, 87.15.Aa, 87.15. Tt 



1 



Overcharging is a phenomenon, happening in many chemical and biological systems. It 
accounts for the situation that more counterions are attracted to the surface of a charged 
macromolecule in an electrolyte solution than needed to neutralize the bare surface charge [1 ] . 
Recent studies show that counterions can form strongly correlated liquid on a charged sur- 
face and the negative chemical potential of the liquid provides an additional attraction of 
the counterions to the surface, leading to overcharging There is a second phenomenon, 
closely related to overcharging, in which the effective charge of a macromolecule inverts 
its sign in an electrolyte solution under certain condition, resulting in the reversal of elec- 
trophoretic mobility sj]. We will call this phenomenon "charge inversion" Q|. In charge 
inversion, all the ions bound to the molecule are considered because they contribute to the 
effective charge, whereas in overcharging, only the charge overcompensation near a molec- 
ular surface is discussed. Without doubt, overcharging occurs while charge inversion takes 
place. However, we do not know for sure if the occurrence of the previous inevitably brings 
in the latter. This is a fundamental question addressed in this study. 

In the past decades, there is a resurgent interest in condensing DNA using charged small 
molecules or salts for diverse applications 5]. This kind of study deals with highly-charged 
polyelectrolytes (PE) and multivalent counterions, which goes beyond the description of the 
classical mean-field theory based on Poisson-Boltzmann equation. Upon addition of multi- 
valent salt, PE is firstly condensed from a solution, due to ionic screening or bridging; the 
condensed PE is then dissolved into the solution when salt concentration is increased to a 
high value f|. Nguyen et al. recently explained these phenomena, called reentrant conden- 
sation, using the idea of overcharging [?J . When the surface charge of PE is reduced to nearly 
zero by condensed multivalent counterions, the correlation induced short range attraction 
dominates the Coulomb repulsion, leading to PE condensation. At a more elevated salt con- 
centration, PE can get overcharged and the Coulomb repulsion dominates; the condensed 
PE thus reenters into the solution. In their theory, overcharging implies the sign inversion of 
the effective charge, and consequently, the reversal of the moving direction of PE in an elec- 
tric field was expected. Reversal of electrophoretic mobility has been observed since decades 
ago 3] and firstly explained by theorists in 1980's [8f\. Nonetheless, the connection with the 
redissolution of PE is still not well understood. Recent study showed that the resolubiliza- 
tion of fd virus bundles happened at high salt concentration without mobility reversa 
Controversy hence arises. Integral equation theories [I], [iq|, II] and simulations 12, 
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have both confirmed the possibility of charge overcompensation next to a PE surface. They 
have also demonstrated that the ion size is an important factor to affect the properties of 
a macromolecule. Moreover, an oscillatory behavior in the integrated charge distribution 
around a molecule has been found in these studies, which indicates the formation of ionic 
layers, alternating the sign of charge. It therefore becomes unclear whether the net charge of 
the whole PE-ion complex inverts its sign or not, even if there is overcharging in the vicinity 
of the surface. Noticeably, theorists also predicted the importance of Bjerrum association 
and the dependence of charge inversion on the ion size 15]. 

The main difficulty in calculation of the net charge comes from no precise information 
about the border of a macroion-ion complex. A simple but accurate way to circumvent this 
difficulty is to do electrophoresis in a weak electric field. In this Letter, we use multiple PEs in 
tetravalent salt solutions as an example and computer simulation as a tool to investigate the 
relationship between overcharging, charge inversion, and redissolution of PE condensation. 
Our system contains four bead-spring chains and many charged spheres, which model anionic 
polyions and monovalent cations, dissociated from the PEs, and tetravalent cations and 
monovalent anions, dissociated from the salts. Each chain consists of 48 monomers and each 
monomer carries a negative unit charge — e. The excluded volume of monomers and ions is 
modeled by a shifted Lennard- Jones potential, U^{r) = £lj [2 (cr/r) 6 — l] 2 , truncated at the 
minimum, with coupling strength e-^j = k^T j 1.2, where k-Q denotes Boltzmann's constant 
and T the absolute temperature. We assume that the monomers, monovalent cations and 
anions have equal effective diameter a = a m and vary the effective diameter of tetravalent 
cation ot from to 4.0a m . This size variation covers a broad range of interest for theorists 
and experimentalists, discussing from point charges to large charged colloids. The virtual 
springs jointing monomers on a chain are described by the FENE potential [lij, £/fene(&) = 
— 0.5/c6^ aa . ln(l — & 2 /fc^j a!E ), with maximum extension b max = 2a m and spring constant k = 
7^Lj/o"m- Solvent is simulated as a medium of uniform dielectric constant e, and hence, the 
Coulomb interaction between two particles of valences z% and Zj is equal to k&T\&ZiZj jr 
where Ab = e 2 / (ineeokBT) is the Bjerrum length and Eo the vacuum permittivity. We set 
Ab = 3o" m . The model is representative of a prototypical PE such as sodium polystyrene 
sulfonate in an aqueous solution at room temperature. The system is placed in a periodic 
cubic box and subject to a uniform external electric field pointing to +x direction. Langevin 
dynamics simulations are employed in this study [17| and the equation of motion reads as 
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mfi = —m^iifi + F[ c ' + ZieEx + rfj(i) where the particles are assumed to have identical mass 
m and rwyi is the friction coefficient. According to the Stokes law, 7, is linearly proportional 



to the particle size af, we thus choose 7, = (a,i/ <j m )T~ l with r = a m Jm/ (k^T). is the 
conservative force acting on i, and fji is the white noise satisfying the fluctuation-dissipation 
theorem. E is the strength of the electric field and we set E = O.Obk-sT/ (ecr m ). The monomer 
concentration is fixed at C m = 0.008<r~ 3 . This C m is well below the overlap threshold and 
the results have been shown representative for the range 0.00032<r~ 3 < C m < 0.016<t~ 3 jl^ . 
In the rest of the text, we use a m , r, e, and k^,T as the units of length, time, charge, and 
energy, respectively. Therefore, concentration will be measured in unit cr~ 3 , strength of 
electric field in unit kBT/(ea m ), and so on. 

We first verified that E = 0.05 is a weak field for our system so that the conformation of 
the chains and the structure of the solution are almost not modified by the electric field. It 
is justified by comparing both the mean square radius of gyration R 2 and the mean center- 
of-mass distance between chains d cm in the presence and in the absence of the electric field, 
shown in Figs. [H(a) and[]Tb). The difference between the two set of data is very small. We 
have verified that the electric field applied here is much weaker than needed to deform the 
chains. 

In Fig. QTa), we observed that the chain size for a t < 3.0 decreases with increasing C s 
up to C*(= 0.002) at which the total tetravalent cations and the 4 PE chains are at charge 
equivalence. A striking phenomenon occurs once C s is increased beyond C*\ R 2 g increases 
or in other words, the chains swell. This chain collapsing and swelling can be regarded as a 
single-chain version of the macroscopic phase separation and redissolution of PEs happening 
in multivalent salt solution ffjj], and have been reported in our previous study 14]. For the 
cases with larger <r t , the collapsing and swelling basically do not take place. It demonstrates 
a strong influence of ion excluded volume on the chain size. In addition to the variation 
of single-chain size, multivalent salt can also provoke subsequently multi-chain aggregation 
and segregation. When the ion size is intermediate (a t = 1.0 and 2.0), d cm curve in Fig. QJb) 
shows a deep valley near C*, with its value smaller than 2R g , which strongly indicates 
chain aggregation. At high C s , d cm increases to the value for neutral polymers; PEs hence 
segregate. It is worth noticing that even there occurs single-chain collapsing in the region 
C s < C*, multichain aggregation may not take place, for example, when o t = 0.0 and 0.5. 

We then studied the integrated charge distribution Q{r) around a PE. Q{r) denotes the 
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FIG. 1: (a) R 2 g and (b) d cm as a function C s for different values of ot indicated in the figures. 
The symbols denote the data obtained in the presence of the electric field, whereas the solid curves 
nearby are the ones obtained in the absence of it. The dotted curve is the data for neutral polymers. 



total charge inside a wormlike tube which is the union of the spheres of radius r centered 
at each monomer center of a PE. For C s < C*, Q(r) increases monotonically from —48 to 
with r, due to electroneutrality JJJ. If C s > C*, a positive peak appears in Q(r) next to the 
PE surface. This peak results mainly from an excess number of the tetravalent counterions 
appearing near the surface; the chains are thus overcharged. The degree of overcharging 
increases with C s . At high C s , Q(r) eventually shows an oscillatory behavior around zero, 
which suggests a repetitive overcompensation of charge inside a wormlike tube. Fig. [2] shows 
an example of Q(r) at high C s = 0.008 for different at- We observed that the positive peak 
near the surface and the oscillatory behavior of Q(r) are both enhanced as a t is increased. 
These phenomena could be explained from the point of view of entropy [18) . Increasing ion 
size decreases the free moving space of ions and thus, decreases the entropy of the solution. 
The correlation of particles therefore becomes stronger, manifested, for instance, in Q(r) 
curve with a more pronounced peak and oscillatory behavior. Notice that for large at, 
monovalent cations can become energetically competitive with tetravalent ones, to condense 
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onto a PE. It is the condensation of monovalent cations which leads to the small bumps at 
r = 1.5 for a t = 3.0 and 4.0. 

The net charge of a PE-ion complex can be calculated by Q e ff = / (ji^ E) , provided 
that the drift velocity in the electric field E and the kinetic mobility //^ of the complex 
are known. Since PE chains move with complexes, we approximated jjtS c ^ by chain mobility 
fi( p \ calculated using Einstein relation: ju^ = / (ksT). is the diffusion coefficient 
obtained by taking time derivative of the mean square displacement of the center of mass 
of a PE chain in y- and indirections and the results are shown in Fig. [3j We found that the 
referenced for neutral polymers decreases with increasing C s due to the jamming effect 
of the salt ions presented in solutions. For PEs, is smaller than the referenced 
while C s < 0.001; it is because extra friction forces act on the condensed counterions, which 
diminishes the chain diffusivity. Since the majority of the condensed ions is monovalent 
in this salt region, the diffusivity is not sensitive to the size of the tetravalent counterions. 
Differently, in the mid-salt region around C*, the main condensed ions are tetravalent; 
depends sensitively on a t . For a t = 0.0, there is no dragging force acting on the tetravalent 
ions; is hence large and the curve shows a hump structure. For at = 1.0, the dragging 
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FIG. 3: I)W as a function of C s for at = 0.0 1.0 and 3.0. The dotted curve denotes the data for 
neutral polymers. 



force takes effect and in addition, chains aggregate and form big complexes, which leads to 
a drastic decrease of the diffusivity. The importance of the latter effect can be seen from 
the absence of the drastic decrease in for o t = 3.0 since the chain aggregation does not 
take place. In the high-salt region, the chains segregate for any a t . The decrease of is a 
result of the combination of the two effects: ionic jamming and PE complexation with large 
number of ions. 

We also approximate by the chain drift velocity . The computed Q e s are presented 
in Fig. HI In the low-salt region C s < 0.001, Q e fr weakly depends on the ion size and 



attains a value rough 
condensation theory 



yequal to —20. This value is smaller than the prediction of Manning 



19] , which states a reduction of the effective line charge density on an 



infinite rodlike chain to — e/A# and yields an effective chain charge —17.2 for our case 
regardless that our PEs are finite and flexible. Q e g turns to become sensitive to a t while 
C s is larger than 0.001. For a t < 2.0, Q e s increases and surpasses zero around C*. It 
clearly demonstrates the occurrence of charge inversion. Q e g attains a value of roughly 10 
at high C s for point tetravalent ions and the value gradually decreases as a t is increased. For 
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<j t = 1.0 and 1.4, a hump appears in Q e s curve and crosses zero line twice, indicating that 
charge inversion takes place only in a window of salt concentration. For a t > 2.0, although 
showing a hump, the curve stays completely in the negative region — no charge inversion 
takes place. Q e g at high C s can be even more negative than in a salt-free solution for large 
ion size such as a t = 3.0 and 4.0. We remark that Q e s obtained here stands for the mean net 
charge of a PE complex which may contain multiple chains, and thus, should not be simply 
interpreted as the net charge of a single PE, specially in the mid-salt region where chains 
aggregate. A general behavior of decreasing of Q c g in the high-salt region was found. The 
onset of the associated phenomenon, the decrease of the electrophoretic mobility at high C s , 
has been observed in experiments [ll|, [21I but never been emphasized before. Both mean 
field theory and integral equation theory fail to predict this decrease. Noticeably, from a 
comparison between Fig. [2] and Fig. HJ we can see that the location of the slipping plane of 
a PE in the notion of electrophoresis depends on the ion size. The smaller the ion size, the 
closer the plane to a PE. The location is roughly at r = [o m + 0"j)/2. 

We remark that the hydrodynamic interaction (HI) was not taken into account in this 
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study. Theorists have shown that HI is screened in a typical condition of electrophoresis 



and thus weak compared to the electrostatic binding force [22]]. This prediction has been 
confirmed by simulations [l^ . In our case, the applied field is weak and the Debye screening 
length is short. Therefore, HI is negligible. Our results show that even if overcharging 
happens next to a PE surface in a high-salt region (Fig. [2]), charge inversion does not 
necessarily come with it. Moreover, a disconnection of the reentrant transition with the 
charge inversion was observed: chains segregate at high C s for at = 2.0 (Fig. Hfb)) but 
charge inversion does not occur. We have specially performed simulations for a t = 2.0 in 
several electric fields of another order of magnitude, ranging from E = 0.005 to 0.1, and 
found consistent electrophoretic mobility and ion distributions around a PE. It confirms 
that the applied field is weak enough so that the following scenario does not happen: PEs 
are stripped of their condensed counterions by the electric field, which leads to the negative 
Q e ff. Even though there is a strong segregation for the case a t = 1.0, the absolute value of 
Q e s is small. These findings suggest that chain redissolution does not come from Coulomb 
repulsion owing to the inversion of the net charge to a non-negligible value. Moreover, 
the reversal of electrophoretic mobility can take place only when the size of multivalent 
counterions is smaU, complying with the theory of Nguyen et a,. Q iu which they assumed 
point ions. Solis and Olvera de la Cruz [15( predicted a different result: no charge inversion 
for small ions. The difference comes from their assumption that the ions have identical size. 
It artificially aggrandizes the effect of ion association when the ions are small. Simulations 
using a model similar to theirs (cf. the second paper in Ref. [jjj]) gave consistent results with 
their prediction. 

We propose the following mechanism for chain decondensation. If C s is high, the increase 
of energy related to the segregation of chains will be small because tetravalent counteri- 
ons present abundantly in the bulk solution and interact with the segregated chains, which 
reduces the energy. Therefore, the entropy increase related to segregation can become a 
dominated term, which lowers the free energy. Consequently, the system favors chain segre- 
gation. The effective charge is another issue. If at is large, the electrostatic correlation will 
be small. The entropy effect related to the ion excluded volume does not produce sufficiently 
strong complexation to lead to charge inversion, although the local overcharging can still be 
possible. 

In summary, we have studied "overcharging" and "charge inversion" of PEs in tetravalent 
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salt solutions. The crucial role of the ion excluded volume on the PE properties has been 
demonstrated. Overcharging does not necessarily bring in charge inversion and the effective 
charge is a nonmonotonic function of salt concentration. Unambiguous evidence has been 
presented that chain redissolution induced by tetravalent salt can take place without charge 
inversion. 
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